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ABSTRACT 

In this paper, we explore instability regions of non-static axial reflection sym¬ 
metric spacetime with anisotropic source in the interior. We impose linear per¬ 
turbation on the Einstein held equations and dynamical equations to establish 
the collapse equation. The effects of different physical factors like energy density 
and anisotropic stresses on the instability regions are studied under Newtonian 
and post-Newtonian limits. We conclude that stiffness parameter has a signifi¬ 
cant role in this analysis while the reflection terms increase instability ranges of 
non-static axial collapse. 

Subject headings: Axial symmetry; Relativistic fluids; Stability 


1. Introduction 

Self-gravitating objects pass through different intense phases of dynamical activities 
during the evolution of the star model. Anisotropy cannot be ignored in the study of rotating 
stars which is closely related with axial symmetry. However, there is no exterior metric 
which coincides with the sources of such interior. It is noted that in conventional celestial 
objects radial and tangential pressures exist instead of purely isotropic fluids. The theoretical 
advances in dicate that such objects in whi ch density ranges upto u < 10 15 gcm~ 3 will be 
anisotropic flRudermanl Il972c ICanutol 1197311 . iBaneriee fe Sanvall (119961 1 classified spatially 
homogeneous axially symmetric spacetimes for an imperfect fluid configuration by keeping 
constant ratio between shear and expansion. 


Herrera et al. (1979J) studied an approach to examine the slow adiabatic contraction of 


anisotropic spheres. They related the radial and tangential pressures by a quadratic law in 
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radial coordinate, defining in this way an anisotropic law (or second equation of state). IRago 
019911 ) explored solutions of the field equations depending upon two arbitrary functions, i.e., 
anisotropic and generating functions which measure the degree of anisotropy and relevant 
physical quant i ties. The solutions are then matched with the Schwarzschild exterior metric. 
Dev fe Gleiserl (1200211 studied the effects of anisotropic pressure on the properties of grav¬ 
itationally bound spherically symmetric object. They found that anisotropic pressure can 
have significant effects on the structure and properties of stellar objects. Particularly, the 
anis otropy can change the critical mass and surface redshift of the equilibrium configurations. 
Mak fe Harkol (2002) presented an exact analytical treatment to the field equations describ¬ 
ing spherically symmetric anisotropic matter configuration. They concluded that anisotropic 
pressure (with radial pressure obeying linear equation of state) increases the maximum radius 
andpnass of the quark star, which in their case is around three solar masses. Dev & Gleiser 
(2003) have extended the formalism developed by Chandrasekhar to discuss the significance 
of anisotropic pressure on the stability of spherical objects against radial perturbations in 
the scenario of Newtonian gravity and general relativity. They have also discussed this for¬ 
malism to study anisotropic sphe res w ith constant energy density and energy densities of 
2 profile. Recently, Mahmoo d et al. 020151) explored spherical collapse and expansion of 
anisotropic cylindrically self-gravitating systems with charged background. 


Chaisi & Maharai (2005) explored a class of exact solutions for anisotropic spherical 


stars with a physically reasonable form of energy density. These solutions hel p to d escribe 
anisotropic nature of compact objects under strong gravitational fields. Hossein et ah (2012) 
discussed the formation of anisotropic compact star with variable cosmological constant and 
ch ecked all the regularity conditions as well as stability of their model. [Sharif _& Bhat ti 


(2014a.b) investigated the role of different physical factors including anisotropic pressure on 
distinct star models. 

The stability of self-gravitating stars is an important issue as only stable equilibrium 
models are viable. A general relativistic treatment is required for a precise evaluation 



axisymmetric spacetime subject to the radial perturbation along sequence of rotating stars. 
Vilenkin & Ford ( 1982 1 investigated gravitational effects for a specific model and showed 


that behavior of the system can be drastically chan ged due to spacetime curvature lead¬ 


ing to stability or instability. Barrow & Ottewill ( 19831) explored the existence and sta- 











































3 


bility of isotropic homogeneous star subject to perturbations in f(R) gravity. iChan et al. 
(119931 . 1994 ) discussed stability analysis via perturbation without using equation of state in 
anisotropic stellar int erior whose results with astrophysical relev ance has also been stud ied 
( Dev fc Gleiser 20031). Recently, stability analysis for spherical ( Sharif fc Yousal 2014a ) as 
well as cylindrical Jsharif fe Yousal 2014b ) configurations are performed by using radial per¬ 
turbation in the framework of f(R) gravity in which they concluded that instability ranges 
depend only on material variables with zero expansion independent of the fact that how much 
the fluid is stiff. However, it shows dependence on the stiffness parameter in the presence of 
expansion scalar. 


Herrera et all (120121) showed that instability range is independent of fluid stiffness with 


zero_ expansion which is compatible with the study of Tolman mass. Babichev & Fabbr 


(120131 ) studied instability of black holes in massive gravity theory and concluded that lin¬ 


ear perturbation around the simplest black hole leads to unstable mode. Roupas (120131 ) 


found a connection between gravity and thermodynamics and discussed stability properties 
under the influence of cosmological constant. We have explored some instability regions 
for self-gra vitating fluids with and wi thout expansion-free c o nditio n by imposing linear per¬ 
turbation ( Sharif fe Bhatti 2014cl ldl). Ixhamesra fe Suneeta ( 2015 ) used the gauge-invariant 


perturbation theory to discuss stability of sphericall y sy mmetric spacetime with anisotropic 
fluid under axial perturbation. Recently, Abbas (2014) discussed dynamical properties of 


commutative black hole and found its total energy. 


Vaik fc Eltgrothl (Il970l ) systematically derived three classes of axially symmetric spa¬ 


tially homogeneous spacetime for an ideal fluid using equation of state. Rao & Neeli ma 


(120131) solved the field equations by using anisotropic features of the universe with axially 
symmetric spacetime and showed that their solution represent evolution of the early uni¬ 
verse. The dynamical stability of rotating or axial reflection symmetric stars against linear 
perturbations as well as the final fate of collapse has not been established definitively. A 
number of interesting attempts i n this direction have already been made with a restricted 
class in a modified gravity (Sharif & Yousaf 2014c, 2015). 


The main idea of this work is to present an analytic treatment in spirit of Chanrasekhar’s 
work (]Chandrasekhar 19641) to identify the instability eras of non-static axial geometry with 
reflection symmetry filled with anisotropic matter. The paper is outlined as follows. I 11 the 
next section, we present some basic equations for axial spacetime which are used to develop 
our analysis. Section 3 provides perturbation scheme to all equations obtained in section 
2 and consequently derives the collapse equation. In section 4, we identify the instability 
regions with N and pN limits from the collapse equation. In the last section, we summarize 
our results. 




























































4 


2. Anisotropic Source and Conservation Laws 

We t ake a n on-static axial spacetime with reflection symmetry in spherical coordinates 
given by (Herrera et al. 2014a,b) 


ds 2 = — A 2 (t , r, 9)dt 2 + B 2 (t, r, 9)(dr 2 + r 2 d9 2 ) + 2 G(t, r, 9)dtd0 + C 2 (t , r, 9)d(j) 2 . (1) 

It excludes explicitly the term representing rotation in the geometry to avoid the complica¬ 
tions in the calculations. We consider anisotropic fluid whose energy-momentum tensor is 
given as 

( 2 ) 


= (fi + P)V a Vp + Pg a p + U a p, 


here 


n a/J = -(II// + 211 j)(K a Kp — -h a p) + — (11/ + 211//) (L a Lp — - h a p ) 


+ n^/^/WLiCi + KpL a ), 


which includes 


Kp = g a p + V a Vp, U KL = K a L?T a p, 11/ = (2A-/C 3 - S a S 0 - L a L 0 )T a p, 
n// = (2L a L d - K a K 0 - S a S 0 )T a p , 

where //. II a/ 3 and h a p are the energy density, anisotropic stress tensor and projection tensor, 
respectively. Also, P is isotropic pressure and II/ ^ II// ^ II KL indicates anisotropic contri¬ 
bution of stress tensor. Moreover, S a , L a , K a and V a are unit four-vectors and four-velocity, 
respectively while a , /3 represent Lorentz indices. In comoving coordinate system, one can 
assume these four-vectors as 


S n = C5i L a = 


\/r 2 A 2 B 2 + G 2 

A 


C k 0 = bsi =-A5» + (3) 


satisfying 


K a L a = S a L a = K a S a = V a S a = V a K a = V°L a = 0, 
K a K a = S a S a = -V a V a = L a L a = 1. 

Using these unit four-vectors, the non-zero components of Eq. (l2jl become 


Too — d ^ 2 , T 02 — ~d^i Tn — ( P + ttH/ ) B 2 , 


T 22 — d 


G 2 (r 2 A 2 B 2 + G 2 


A 2 V' 


A 2 


P + 3 n „ 










-5 - 


B 


Tv2 = TIkl [ ^Vr 2 A 2 B 2 + G 2 ) , T 33 = 


P ~ 3 (n/ + Tin) 


C 2 . 


( 4 ) 


In order to describe dynamical nature of any self-gravitating system, the conservation 
law, T a ^.g = 0, has a crucial role which from Eqs.([2]) and (J4]) for a = 0,1 yield the following 
couple of equations 




x 


B | G 1 

B + G + r 2 A 2 B 2 + G 2 


r 2 AAB 2 + GG + r 2 A 2 BB 


A 2 B 2 


x 


r 2 A 2 B 2 + G 2 
B C 


2 B C 


G 2 


B + Cl + A 2 B 2 \ B ' G A 


B G A C 
+ ^~^ + G 


+ (n + P) 


B C 


+ 


n 


a 


r 2 A 2 B 2 + G 2 


r 2 A 2 B 2 


B C 
B~C 


+ G 2 


G 

A 

c) 

U‘ 

A 

C 


= 0, 


p' + - (2n' + u' n ) + 

'A! 2 B’ 2 

X 1 ~A + ~B + r~ 


-P + — (211/ + II//) 

y 


r 2 A 2 B 2 


C’ 3GG' 

~C + 2 + r 2 A 2 B 2 + G 2 


{rB)' 

rB 


r 2 AB 5 


6P, G, 4GG 0 

X ' A + B + C + r 2 A 2 B 2 + G 2 ' r 2 A 2 B 2 + G 2 \ A 


(■ r 2 A 2 B 2 + G 2 ^ 
Ar 2 A 2 B 2 

+ 


-TIkl, e 


r 2 AB 5 


(r 2 A 2 B 2 + G 2 )% 

^ + §Y>im 


+ 


fir A A A B 4 


(r 2 A 2 B 2 + G 2 ) 2 
= 0, 


55 + ^ _ ^ ^ _ _p r 2 A 2 G 2 B 2 


G’ (rB)’ 

A r 2 AB 2 ) ( r 2 A 2 B 2 + G 2 ) 2 V 2G + rB 


( 5 ) 


( 6 ) 


where dot and prime stand for differentiation with respect to t and r while subscript 9 
indicates ^-differentiation. 


3. Perturbation Approach 


Here we use perturbation technique to obtain perturbed form of all the previous equa¬ 
tions up to first order keeping the perturbation parameter e in the interval (0,1). The initial 
configuration of the system is considered to be static while after perturbation it enters into 
non-static phase with the same time dependence of metric coefficients. Consequently, the 


metric and matter variables are perturbed as follows (IHerrera et al.l 120121: ISharif fe Yousal 


2 014 a.b) 


A(t,r,9) = A 0 (r,9) + eT(t)a(r,6) 


( 7 ) 









































6 


B(t,r, 9) 

= B 0 (r, 9) + eT(t)b(r, 9), 

(8) 

C(t, r , 9) 

= C 0 (r,9) +eT(t)c(r,9), 

(9) 

G(t, r , 9) 

= G 0 (r,9) + eT(t)g(r,9), 

(10) 

At, r, 9) 

= Ho (r,6)+eji(t,r,6), 

(11) 

P(t, r, 9) 

= P 0 (r,9)+eP(t,r,9), 

(12) 

n j(t,r,0) 

= fl/ 0 (r, 9) + £lT/(f, r,9), 

(13) 

n u{t,r, 9) 

= II//o(r, 9) + dd//(t,r, 9), 

(14) 

n KL (t,r,0) 

= B. KL0 (r,6) + AlKL(t,r,0). 

(15) 


Using the above equations, the first conservation law ([5]) is perturbed as 


h = 


ho ^ TV" + 7T + V (r 2 aA )£5 + hG 0 + r ~bB 0 Afy }> + (710 + P 0 ) 
-do n-o Zq 



x ^i<!r 2 
6 

X ' A) 


26 2c 


, G l (Jl + A _ A + JL 

AlBl \ B n Go A 0 Co 


c 

C~o 


+ 


n 
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3Cfi 


r 2 AlB% 


0-^0 \ J -'0 Oro 
b C 


n 


70 


B 0 G 


+ c: 


_g_ 

Gq 


a 

Ao 


c 

Co 


T. 


(16) 


It is interesting to mention here that only the non-static part of first conservation law exists 
while the static part vanishes for this case. Similarly, the static part of second conservation 
law leads to 


P'o + g ( 2n /o + n //o) + 


Pq + ~ (2II/o + n //0 ) 

y 


Cp 3G 0 Gj rM IBS 
C„ 2 Z» 


X 


A' 0 2 B\ 


— + 

A-o Bq 


+ 2_I_^o 
r r B n 


r 2 A 0 Bl 

3 

Z% 


n 


r 2 A 0 Bl (A 


KLOfi 
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i06> 


An 


6 B, 


09 


Bo + Co + 


C oe , 4G O G O0 , Ar 2 A 2 Bl {A 


+ 




Ao 


+ 


B, 


09 


Bo 


n 


KLO 


fioAA^Bl fA' 0 GqAqo \ ^ 0 r 2 A 2 G 2 B 2 


Z 2 o 


Ao CAoBl) 


zl 


G’o IB’ 

-2- h -h — = 0. 

2G 0 r + B 0 ' 


( 17 ) 


The non-static part of the second conservation law after perturbation turns out to be 


B 2 
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-Ap' + - (211; + n' 7J ) + --AP + -(211/ + n n ) ^ + 


Bl 


2, - 


C' 3G 0 G' 


Co 


aa^Bq ( a'q ±j 0 
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x 


x 


Ape 6-Bpg C 0 q 

AT “Bo" AT 


4G 0 G 0 , , 4r 2 A 2 S 2 /A o0 


+ 




An 


An 
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0^00 
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B' 
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T 
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Bp 

2b 
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L -°0 
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p.r 4 A^ 




A + 


c 
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r r Bp 
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3 


n 


ifLO,0 


a 36 3Z 
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Ao Bp 

r 3 ApBqTIklp 
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Ao 
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3 b 
B( 

r 
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3 Z\ f Ap e 6 Bpe Cpe 4G 0 G 0 e 4r 2 AyBy f A oe 

Aq Bp Cp Zp Zp 
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3 A 0 Sq 
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Zj 


6 Bpe 
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Bpe Bp 


a c 

Ap Cc 
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Ao 
4G O G O 0 


0 / 0 


x 


9 9e 
Gp Gpe 


9o r4 Ap ( f o_ 

Zl lUo 


Z 0 + 


4r 2 A 2 S 2 (2a . 2b 


Zp 

Go Ape 




r 2 A 0 B 2 \ G, 


\ Ao Bp Zp 

ae a 

Ap 


a 


Ap Br 


4- + 


A 


00 


jipGlAlr 2 f2g_ 2a _ 2Z 

Z 2 V Go Ao Z n 


+ 


AioGgA^r 2 ( # 


7 2 

z o 


+ 


26 \ 

BoJ 

B~p 


0/ 0j 


= 0 . 


(18) 


The non-static part of 02-component of the field equation as obtained in Eg. (I A3 1) can be 
written as 

IT + mT + nT = 0, 

here l, m and n are functions of r and 9. The solutions of this equation involve stable 
and unstable configurations. To describe t he instability range , we concentrate over t he 
unstable part which h as the following form (IHerrera et al.l 120121: ISharif fe Yousafl l2014ailb ; 
Sharif &; Bhatti 2014allbl ) 


T(t) = — exp (\/at), where a = 


—m + \J m 2 — 4 In 

21 ' 


(19) 


For the solution to be real, we take a > 0 with certain constraints that m < 0 while ( n > 0. 
Such a solution describes static system which undergoes collapse with large past time. 

























































Now, we calculate non-static part of anisotropic stresses in terms of static matter 
profil es to evalu ate t he instability regimes. We take the equation of state introduced by 


Harrison et all (1 1 9650 which relates pressure with energy density using adiabatic index Y as 
follows 

P = r-^y a, (20) 

ho + M) 

where Y represents rigidity or stiffness in the fluid which is taken to be constant in our 
stability analysis. The value of ft from the non-static part of first conservation law after 
perturbation can be obtained by integrating Eq. (fT6l) with respect to t as follows 


h = - 


ho 


— + — + — ( r 2 aA 0 By + gG 0 + r~bB Q Afy + (ho + P 0 ) 
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Using this value of g in Eq. (l20lh we obtain 

^ n 70 
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( 21 ) 


here 
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Substituting these values in Eq. (fT8]h it follows that 
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where 


Z 0 = r 2 A 2 0 B 2 + G 2 , 


Z = 2 r 2 A 2 0 B 2 


a 6 
Aq Bq 


+ 2gGo- 


(23) 


This is the required collapse equation with the constraints Pg, n 70 , P' II0 < 0 which is very 
useful to investigate the instability regions for our systematic analysis. 


4. Vorticity Tensor and Instability Regions 

This section investigates vorticity tensor and dynamical instability ranges for non-static 
axial spacetime with the help of equations obtained in the previous section particularly 
the collapse equation subject to N and pN limits. The role of stiffness parameter and its 
dependence on physical factors are also analyzed in this scenario. 
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The kinematical variable responsible for producing local spinning action of anisotropic 
fluid configurations is the vorticity tensor. For reflection axisymmetric spacetime, this tensor 
in terms of four vectors, K a and L a , can be expressed as 


where 


n 


G (G' 2A'\ 

2 By/Z U A ) ' 


For the vanishing of vorticity scalar, either G — 0 or ^ — BA = 0. If we take G = 0 then it 
leads to the vanishing of vorticity scalar. On the other hand, if we take c -jj — = 0, then it 

vanishes the metric coefficient A{t , r, 6) describing the temporal component of the spacetime 
as follows 


G' _ 2 A! 
~G~^~ 



= 0, 


= 0, 


where C = C(t,6) is an arbitrary function of integration. Consequently, GC = A 2 , which 
implies that for G = 0, we have A = 0 disturbing the existence of our non-static axial 
spacetime. Hence, 


G' 

~G 



Thus, we take G = 0 with regularity condition at the center indicating that vorticity of 
axisymmetric spacetime exists if and only if its reflection degrees of freedom exist or more 
precisely G = 0 G = 0. Consequently the assumption G = 0 in the dynamical evolution 
of non-static axisymmetric anisotropic metric gives zero value to non-diagonal scale factor, 
G. whose dynamics has already been discussed in GR (ISharif fe Bhattil 12014bh as well as in 
modified gravity theory (ISharif fc Yousall 2014ch 


4.1. Newtonian Limit 

For N limit, we take Ho = 1, Bq = 1, Co = r, Go = r so that Zo turns out to be r 2 for the 
instability analysis. We also discard the terms of order where nio is the static profile of 
the mass function. The physical requirement of the collapsing matter, i.e., Pq, H/ 0 , n// 0 < 0, 
is also imposed in N approximation. By making use of the above mentioned constraints, the 
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collapse equation yields 
TT 
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The system will be unstable until it satisfies the following relation 

f (26' + ? + £ - o') + & + Ai 


r < 


„/ 1 11 „ 1 n ^klo, 9 (ou 1 2e'\ 

V + TrV + Si - -3v 5T ( 3f> + vj 


(24) 


0 


where 


V= Ub + — + -'j &, 

\ r r J 

while remaining quantities are defined in Appendix A. It is well-known that instability will 
emerge as long as all the terms given in the above inequality are positive. For this purpose, 
we need to take |Al|, \Gi\ and |^| instead of Ai, Gi and Q?. We find that the adiabatic 
index depends on the static profile of matter variables of axial geometry. 


In the above instability constraint, the quantities \Gi\ and |£? 2 | incorporate meridional 
effects that arise due to non-zero vorticity vect or of the collapsing system. It is well- 
known from the work of Herrera et ah ( 2014ah that invoking of reflection effects in 
axially symmetric anisotropic stellar object causes the emission of gravitational radia¬ 
tions. These radiations induce the loss of both energy and angular momentum, which 
consequently boosts up the instability of the reflectional axisymmetric body. 

The quantity A\ includes anisotropic contribution of axial geometry. It is seen from 
the expression (1241) that anisotropy tends to produce complications in understanding 
its role in the stability of axial systems. However, if one considers positivity of all terms 
in denominator and numerator of the above expression then it is seen from (12TD that 
anis otrop ic press ure tends to increase instability regions. This result is well-consistent 
with Chan et ah (1993). 
























4.2. Post-Newtonian Limit 


For the instability era in the pN approximation, we assume 


A 0 — 1 


m 0 

) 

n 


5o = l + — , 
r i 


(25) 


and the terms of the order — while discarding the terms containing higher orders of —. 
The system will be unstable in the pN region if it satisfies the inequality 


r faC + Gs + Ai + X, 

A:j 


(26) 


where static profile terms X\ and Q 3 are non-diagonal and diagonal components of scale 
factors at pN epoch, respectively, while A 2 and A 3 incorporate anisotropic effects in the 
evolutionary phases of collapsing self-gravitating axial stellar object. These terms are given 
in Appendix A. 


4-2.1. Restricted Class of Anisotropic Axial Spacetime 


On assuming Q 3 = 0, our instability constraint at pN approximation of axisymmetric 
object with reflection symmetry reduces to 


r < 


IfCC + A2 + Xi 


A3 


(27) 


This describes instability range of the restricted class of non-static axial geometry since it ex¬ 
cludes explicitly rotations around the symmetry axis, i.e., dtdcj) a s we ll as the refle ction terms. 
This result coincides and supports already calculated solution (iSha rif fe Bhatti 2014b! ). 


4-2.2. Reflection and Restricted Class of Isotropic Axial Spacetime 


On taking equal all principal stresses as well as zero value to Q 3 , one can find insta¬ 
bility regions of restricted class of isotropic axisymmetric spacetime from expression (J26]). 
However, apart from that by assuming only first of above limits, one can get dynamical 
instability constraint of reflection axisymmetric compatible with perfect fluid. All possible 
stellar models of reflection axial symmetric system coupled with perfect (isotropic) matter 
configurations have have been explored in detail by Herrera et ah) (12015(1 . 


We see that the adiabatic index T plays a central role to investigate dynamical instability 
of the relativistic system. It is worth mentioning that for T < | and T < 1, the spherical and 
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cylindrical relativistic objects become unstable respectively thereby enforces the importance 
of index T. Infact, the adiabatic index also known as stiffness parameter demonstrates 
how much relativistic fluid is stiff. We have established the relevance of such index in the 
dynamical instability as seen from expressions (12T11 and (12()jl depending upon the static profile 
of the structural properties of the system. The system would be in complete hydrostatic 
equilibrium, if (during evolution) adiabatic index is able to attain value equal to the right 
hand side of the expressions given in (T24|) and (l26lh However, if stiffness parameter attains 
a value greater than the right hand side of expressions (1241) and (1241) . then the relativistic 
system begins to move in the stable window, thereby ceasing the collapsing mechanism. 


5. Conclusions 


This paper is devoted to investigate dynamical instability of non-static axially symmetric 
spacetime by choosing reflection term in the geometry. Since rotating stars are more stable 
than non-rotating, so for the instability regions, we have neglected the term representing 
rotation in the general non-static axial spacetime. It is worth mentioning that for axially 
symmetric sources perfect fluid distribution seems to be inflexible restriction, even in the 
static case. On the other hand, Bondi coordinates are known to be very useful for the 
treatment of gravitational radiation in vacuum, but are not particularly suitable within the 
source. An analytical approach, which shares some similarities with ours, although restricted 
to the perfect fluid case, can be found in the literature. Therefore, here, we have considered a 
source which includes all nonvanishing stresses compatible with the symmetry of the problem 
to carry out our systematic analysis. 


We have explored the field equations and corresponding conservation laws in this sce¬ 
nario. We have found three independent components from the conservation l aw w hile there 


exist onl y two co mponents in the case of spherical and cylindrical spacetimes (IHerrera et al. 


2012 ; Sharif fc Bhatti 2014a b 30)- The radial perturbation is used for metric as well as 
material variables to obtain perturbed form of these dynamical equations. We have ex¬ 
plored static and non-static parts of independent components of the conservation law. It is 
found that only non-static part for the first conservation law exists and static part vanishes 
while the remaining equations have both static as well as non-static components. Using 
02-conrponent of the field equations with perturbation technique, we have found a solution 
which corresponds to both stable and unstable configurations and star t collapsing at large 
past time diminishing its areal radius ( Sharif fe Yousaf 2014a b 3, 2015 b 


We have developed a general collapse equation to examine the instability regions using 
non-static parts of anisotropic stresses and the solution (fT2jh We have explored two insta- 
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bility ranges under N as w ell as pN limits and found that instability range is defined by the 
adia batic index (ISharif fc Yousafl l2014d. 120151) unlike expansion-free case [where it has no 


role ([Herrera et al. 20lJ)]. The adiabatic index depends upon static profile of the energy 


density, anisotropic pressure and the reflection term in the spacetime. We conclude that 
reflection symmetry increases the unstable range of the axial geometry. The system will 
remain unstable until it satisfies the relations (|24|) and (1261) while their violation will lead to 
stable configuration of the model, ft would be interesting to examine the role of dissipative 
terms like heat flux on the stability of non-static axial geometry. 


Appendix A 

The 02-component of the Einstein tensor corresponding to our line element in Eq.([T]) 
takes the form 
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For a = 3, the conservation law, T ' ^ = 0, leads to the following equation 
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Using Eqs.(j3}, (IA1D and (JTj) - (fTTi) . the non-static part of the 02-component of the held equa¬ 
tions, G a /3 = 87rT a i3, becomes 
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The quantities introduced in the unstable range in pN limit are given as 
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